Trapping Horizons as inner boundary conditions for black hole spacetimes 



am6n, 3 '@ and Jose Maria Ibanez 3 'H 



o 
o 

(N 

in 

00 



o : 

cr 

i , 

5h . 
W), 



> 
(N 

00 
(N 

o\ 

O 

o 



X 



Jose Luis Jarairiillo, 1,2: B Eric Gourgoulhon, 2 '0 Isabel Cordero-C 

1 Institute) de Astrofisica de Andalucia, CSIC, Apartado Postal 3004, Granada 18080, Spain 
2 Laboratoire Univers et Theories (LUTH), Observatoire de Paris, CNRS, 
Universite Paris Diderot ; 5 place Jules Janssen, 92190 Meudon, France 



Departamento de Astronomia y Astrofisica, Universidad de Valencia, 

(Dated: 18 September 2007) 



Valencia, Spain 



We present a set of inner boundary conditions for the numerical construction of dynamical black 
hole space-times, when employing a 3+1 constrained evolution scheme and an excision technique. 
These inner boundary conditions are heuristically motivated by the dynamical trapping horizon 
framework and are enforced in an elliptic subsystem of the full Einstein equation. In the station- 
ary limit they reduce to existing isolated horizon boundary conditions. A characteristic analysis 
completes the discussion of inner boundary conditions for the radiative modes. 
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General problem. The aim of this report is to discuss a 
set of inner boundary conditions (BC) for dynamical evo- 
lutions of black hole spacetimes using an excision tech- 
nique. These BCs are derived in the context of the dy- 
namical trapping horizon framework 0,0, 3, 4]. In paral- 
lel with the recent black hole numerical studies based on 
free evolution schemes, which have led to the successful 
simulations of binary black hole coalescence through the 
merger phase (see e.g. |5[ for extensive references), a 3+1 
scheme for a fully-constrained evolution of Einstein equa- 
tion has been presented in Ref . Q . This approach maxi- 
mizes the number of elliptic equations to be solved during 
the evolution resulting in a coupled elliptic-hyperbolic 
PDE system pj ■ Spectral methods 0] are then employed 
both to solve the elliptic subsystem and to handle the 
spatial part of the relevant hyperbolic operators. We deal 
with the black hole singularity by means of the excision 
technique. This raises the question about the appropriate 
choice of inner BCs on the excised sphere, both for the 
elliptic and the hyperbolic parts of the system. Regard- 
ing the hyperbolic equations, this inner boundary issue 
is intimally related to the metric type of the world-tube 
hypersurface generated by the time evolution of the ex- 
cision sphere. As observed in Ref. certain choices 
for the excision surface render this excision hypersurface 
partially time-like, leading to ill-posedness if inconsistent 
BCs are supplied for the radiative modes. A solution 
to this problem is suggested by the quasi-local approach 
to the evolution of black hole horizons, embodied in the 
dynamical trapping horizon framework (see review arti- 
cles 0, HJ and also Ref. [HI). This formalism moti- 
vates a natural geometric choice for the excision surface. 
The basic underlying idea goes back to Eardley's work 
and consists in modeling the black hole horizons by 
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. world-tubes sliced by apparent horizons, that sat- 



isfy certain additional conditions guaranteeing the phys- 
ical growth of the horizon area (see below) . On the one 
hand, apparent horizons at each given 3-slice of the time 
evolution provide non-ambiguous geometric choices for 
the excision sphere that are guaranteed to lay inside the 
event horizon, and therefore are causally disconnected 
from the rest of the spacetime. On the other hand, dy- 
namical trapping horizons are space-like hypersurfaces 
suggesting that no conditions must be supplied at the 
inner boundaries for the modes propagating in the bulk. 
In sum, this proposal recasts Eardley's program [ll[ in 
the dynamical trapping horizon setting. In the following 
we describe the fully-contrained scheme, then we present 
inner BCs for the elliptic part that guarantee that the 
excised sphere generates a (dynamical) trapping horizon, 
and finally we show that the combination of a Dirac-like 
gauge Q and dynamical trapping horizon inner BCs for 
the elliptic part of the PDE system, actually imply that 
no BCs must be prescribed for the hyperbolic fields at 
the inner excised sphere. 

Fully- constrained evolution scheme. In the setting of the 
standard 3+1 decomposition of a spacetime (M,g) by 
spatial slices (St), Ref. @| proposes a particular initial- 
boundary problem for the spacetime evolution from an 
initial Cauchy slice. Let us denote by n the unit time- 
like normal vector to S t , the spatial 3-metric by 7, i.e. 
7 = g + n <g> n, and define the extrinsic curvature of S t 
as K — —T;C n ~f. The evolution vector t = d t is de- 
composed in terms of the lapse function N and the shift 
vector f3, as t = Nn + f3. In addition, we introduce a 
fidutial flat metric /, satisfying Ct f = dtfij = 0. Now 
we proceed by performing a conformal decomposition of 
the 3+1 fields: 7 = * 4 7, K = ^ 4 A + |i^7, where 
K = 7 y 2£y, the representative 7 of the conformal class 
of the 3-metric is chosen to be unimodular, i.e. satisfies 
det(7) = det(/), and the traceless part A of K is writ- 
ten as A ij = ^ (t) l fP + - \b k {3 k ^ + <9 f 7 y ) , D 
being the Levi-Civita connection associated with 7. In a 
second step, a coordinate choice must be adopted. Fol- 
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lowing the prescriptions in Jg|, namely maximal slicing 
and Dirac gauge, we set 



(1) 



where V is the Levi-Civita connection associated with 
the flat metric / (see Ref. [l2| for a discussion and re- 
lation to other coordinate choices). Conditions (Q]) fix 
the coordinates up to boundary terms. The Dirac gauge 
condition will play a key role in the following, whereas 
maximal slicing can be relaxed to an arbitrary K vanish- 
ing asymptotically near spacelike infinity. Inserting the 
conformal decomposition and gauges ([1]) into Einstein 
equation results in a coupled elliptic-hyperbolic system 
|6|. The elliptic part can be written as 
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D k D k i + ■=D i D k p k + 3 ff k /3 k = 
D k D k N + 2£) fc ln* D k N = 



S9[%N,0,t\ 

Si[*,N,0,j} (2) 
S N [N,%0,i] , 



where the first equation on W follows from the Hamilto- 
nian constraint, and the equation for the shift f3 results 
from the simultaneous imposition of the preservation of 
the Dirac gauge in time, i.e. dt(T> k j kl )=0, together with 
the momentum constraint. The Dirac gauge ensures the 
elliptic character of this equation. Finally the third equa- 
tion follows from dfK = 0. Stp, Sp and Sn represent non- 
linear sources given in Ref. @. Note the similarity with 
the extended conformal thin sandwich elliptic system fl3l ] 
for the construction of initial data. In the present con- 
text, Eqs. @ are meant to be solved along the whole 
evolution, not only on an initial slice. Regarding the evo- 
lution part, we solve for the deviation h of the conformal 
metric from the flat fidutial one /, i.e. h = j — f. We 
choose a second-order form for the evolution equations, 
that can be formally written as 



-mt-W 1 VkVlh 
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■LeLpbP = S% (3) 



where the nonlinear sources SV [N, ty, j3, 7] do not con- 
tain second derivatives of h. Eqs. ((2|) and ([3]) are solved 
in Ref. [|| inside a spacetime region bounded by an outer 
timelike tube at large spatial distances. We focus here 
only on the inner BC problem. On a first stage, dy- 
namical trapping horizon considerations will provide in- 
ner BCs for the conformal factor 'J, the shift (3 and the 
lapse N . In a second step we will analyse the hyperbol- 
icity of the subsystem and, most importantly in the 
present context, we will evaluate its characteristics fields 
and speeds to assess if inner BCs must be provided at all 
for h. 

As mentioned above, we do not discuss here the impor- 
tant outer BC problem. In this sense, a very interesting 
alternative has been recently presented by Moncrief et 
al. They propose a (conformal) 3+1 constrained 

scheme, which differs crucially from 6] in one point: the 



chosen slicing, involving constant mean curvature slices, 
extends up to future null infinity X + , a natural boundary 
for physical outgoing radiation conditions. This strategy 
permits to bypass the boundary problem at the outer 
timelike border. The feature of [lj] we highlight in the 
context of the present work is the shared adoption of 
an inner excision approach to the black hole singularity 
problem. An alternative geometric choice for the inner 
surface is proposed in [lj|, namely the use of minimal 
surfaces. However, our proposal of rather employing ap- 
parent horizons instead, straightforwardly translates also 
into their scheme. 

Inner BCs for the elliptic part: dynamical trapping hori- 
zons. Quasi-local approaches to black hole horizons aim 
at modeling the boundary of a black hole region as world- 
tubes of apparent horizons (St). At each point of a given 
spacelike closed surface St we can define (up to total 
rescaling) two null vectors t and fc, satisfying k ■ £ = —1 
and spanning the plane normal to St- Denoting by q the 
metric on St induced by the ambient metric g and by 
es the associated area element, we can define the expan- 
sion 0( v ) and shear cr^ along any vector v normal to 
S t by C v e s = 9^e s and = C v q- 9^q. The 

surface St is trapped (l5j if light rays emitted from it 
locally converge: 6>( fe ) < and 6& < 0. In the lim- 
iting case in which one of the expansions vanishes, St 
is called a marginally trapped surface (MTS). Since we 
will deal with asymptotically flat 3-slices, we can unam- 
biguosly define an outgoing null normal, say £, as the 
one pointing towards spacelike infinity. Then, condi- 
tion 9^ — defines a marginally outer trapped sur- 
face (MOTS) 0. In contrast with MTSs, MOTSs im- 
pose nothing on 9^ k K Apparent horizons are outermost 
MOTSs. In this context, quasi-local dynamical trapping 
horizons 7i are S 2 xl hypersurfaces sliced by MOTSs 
(St) and satisfying 9^ < 0. Actually, slices (<S t ) are in- 
deed MTSs but, motivated by inner BCs below, we wish 
to stress the underlying MOTS structure. Following Hay- 
ward [l|, 7i is & future outer trapping horizon (FOTH) if, 
in addition, Ck 6^' < holds. This represents a stabil- 
ity locally outermost condition, essentially stating that 
the interior of H is a trapped region. FOTHs can be 
either null o spacelike hypersurfaces, the former repre- 
senting stationary situations and the latter dynamical 
ones. Alternatively, dynamical horizons (DH) introduced 
by Ashtekar and Krishnan [2j substitute the condition on 
Ck 9^ by the requirement of 7i to be spacelike, station- 
arity being represented by (null) isolated horizons (IH). 
Both in FOTHs and DHs, condition 9^ < guarantees 
that the horizon area is never decreasing. In the dynam- 
ical contextjFOTHs and DHs have been shown to be 
equivalent [l(| G3]- 111 our 3+1 description, slices (St) 
of TL will always lay within a spatial surface St of the 
chosen 3 + 1 slicing. Denoting by s the unit spacelike 
normal vector to St laying in S t and pointing towards 
spacelike infinity, we can perform a 2+1 decomposition 
on the horizon. In particular, the metric q induced on St 
can be written as q = 7— s®s and the shift can be decom- 
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posed in its normal and tangential part as: (3 = ft s — V, 
with /3 X = /3 • s and V • s = 0. 

A most important result in this context is the foliation 
uniqueness theorem by Ashtekar and Galloway [18| stat- 
ing that, for a given DH Ti, there exists a unique folia- 
tion (St) by MTS's. Using this, we can define a canoni- 
cal vector h as the vector tangent to Ti, normal to each 
St and that Lie-drags each MTS St of Ti into another 
one S t +st- It constitutes a natural evolution vector on 
Ti and can be decomposed as h = Nn + bs, where the 
normalization N follows from requiring St € S t and is 
fixed up to a factor only depending on t. Defining a pa- 
rameter C as (half) the square norm of h with respect 
to g, i.e. C := h ■ h/2 = b 2 — N 2 , it follows from the 
above-commented metric type of FOTH's that C > 0; 
strict inequality b — N > holds in the DH situation and 
b — N — in the equilibrium (null) IH case; accordingly, 
we normalize the null vector £ as the limit of h in the 
stationary case: £ = N(n + s) [ToL [20l| - 
Our criteria for setting BCs for Eqs. @ are: a) to en- 
force the excision world-tube Ti to be sliced by MOTS, 
and b) to recover IH BCs [H, ES EH E3, El at the equi- 
librium limit C = 0. Motivated by this second point, 
but ultimately justified by the inner boundary analysis 
of Eqs. we choose a coordinate system adapted to Ti 
by demanding t to be tangent to Ti. This implies ft 1 - = b, 
and we have 



h = t + V , ft 1 - - N>0 



(4) 



i) Geometric conditions for Ti. The first two BCs are 
provided by 1) the geometric definition of St as a MOTS: 
6>W = 0, and 2) the Lie-dragging of MOTS into MOTS 
by h inside Ti (trapping horizon condition): Ch 0™' = 0. 
The first one yields 



4s • Din* + t> s + ^~ 2 K(s, s) 



^ 2 K 







(5) 



where tildes refer to the conformal metric 7; in partic- 
ular, s — ^> 2 s. The second geometric condition follows 
from the projection onto St of one component of Einstein 
equation and results in the elliptic equation [Tl| 

[- 2 A - 2L ■ 2 D + A] (ft 1 - — N) = B(/3 ± + N) , (6) 



where L 
B 
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K kl s k q l t ,A= ± 2 R- 2 D-L-L-L-8TrT(i,k), 
4irT (£,£), T is the stress-energy ten- 
sor, £ = n + s,k = (n- s)/2, and 2 D, 2 A and 2 R 
are respectively the covariant derivative, Laplacian and 
Ricci scalar of (St,q). The non-negative character of 
the rhs term in ©, together with the lhs elliptic op- 
erator under the FOTH condition (closely related to the 
stability condition in [13)) guarantees the positivity of 
(ft 1 - - N) in Moreover, null- like condition ft 1 - = N 
[ljJE3] is recovered in the stationary IH limit, for which 
a-W — = T(£,£). Condition © provides a relation be- 
tween combinations (ft 1 - — N) and (ft 1 - + N) : given one, 
the other is fully determined. 



ii) Gauge conditions for the tangential part of the shift. 
Let us express the shear tensor along h, cr^ h \ using the 
coordinate system adapted to Ti: 



2a 



(h) 



dqij d 



( 2 D i V j + 2 D j V i - 2 D k V k q ij ) 



(7) 



Imposing as a coordinate choice the vanishing of the first 
parenthesis in the rhs results in 



2 AVj + 2 DjVi - 2 D k V k qij = 2a. 



(h) 



(8) 



an elliptic equation whose source is determined by the 
evolution equation of the shear cr^ (tidal equation): 



jCho-W - -q\q l / m i n w m knl ~c 2 q k l q l J k m k n w m knl 



8irC 



-X<f V T^)qii 



, (9) 



where W is the Weyl tensor. Condition ([5]) fixes the 
tangential part of the shift V up to a linear combination 
the six conformal symetries in the kernel of the elliptic 
operator in the lhs. We determine this conformal sym- 
metry in the evolution by continuity with the conformal 
Killing symmetry prescribed on the initial data. In the 
stationary limit, where h tends to £ and er^ = 0, the 
vanishing of the rhs in Eq. ([5]) leads to the conformal 
Killing condition on V and, given the rescaling proper- 
ties of the conformal Killing operator, IH condition for 
V in Refs. [H, ES El is recovered. 
Hi ) Slicing condition. Combined results in Refs. [TtI 
Il8j show that, for different choices of 3-foliation (S t ), a 
given MTS St on a given initial 3-slice evolves generically 
into distinct DHs. However, all these DHs are ultimately 
expected to approach the event horizon, and therefore 
there is no preferred candidate on the sole basis of the 
dynamical trapping horizon framework. The choice of 
inner BC for AT, must be adopted on the basis of the 
well-posedness of the elliptic-hyperbolic system and the 
specific numerical needs. In practice, this issue must be 
numerically addressed. Having said this, Eq. ([6]) suggests 
an alternative in this context: prescribing an inner BC 
for (ft 1 - — N) determines (ft 1 - + N) algebraically. Such is 
the case of the proposal in [24| , where the choice of that 
DH locally maximizing the area rate of change of the 
slice S t leads to: ft ± — N = —const • 9^ k \ with const > 0. 
Note that if, alternatively, inner conditions are provided 
for (ft 1 - + N) [resp. N], then Eq. ^ must be solved as 
an elliptic equation on St for (ft 1 - — N) [resp. ft- 1 ]. 
Inner BCs for the hyperbolic part. Assessing the free- 
dom in prescribing inner BCs for Eqs. ([3]) is a key step 
in the implementation of the fully-constrained evolution 
scheme. A first analysis of the general issues concerning 
hyperbolicity in Eqs. ([3]), has been carried out in Rcf. 7\ 
by writing down the evolution equations as a first-order 
system in conservative form, i.e. dtU + A l (U)diU = 
F[U,...], where the evolving variable vector U is given 
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by U = (h,d t h,T)h) and matrices A 1 are straightfor- 
wardly determined from Eqs. ([3]). First, it is shown that 
imposing the Dirac gauge in {T]) indeed guarantees the 
real character of the eigenvalues corresponding to ma- 
trices A 1 , and therefore the hyperbolicity of the evolu- 
tion system. Of particular relevance for the present in- 
ner BC discussion is the explicit determination of the 
(non- vanishing) characteristic speeds associated with the 
vector s normal to the excision surface St, resulting in 



= -/3 1 - ± N (each one of multiplicity 6). (10) 

Taking into account the inequality in ((¥]), consequence 
of the choice of a coordinate system adapted to the DH 
TL by enforcing condition ([6]), we conclude the absence 
of ingoing radiative modes into the integration domain 
E t at the excision surface. Therefore no inner BC what- 
soever must be prescribed for the hyperbolic part, as a 
consequence of our choice of BCs for the elliptic part. 
This confirms our initial motivation for using space-like 
excision worldtubes in the evolution and shows the key 
interplay between elliptic and hyperbolic modes in the 
coupled fully-constrained evolution system. 
Discussion. In the context of constrained schemes for 
excised black hole evolutions such as Rcfs. inner 
BCs ([5]) and ((6j> , together with the essentially free choice 
of 3-slicing, characterize the inner excision hypersurface 
Ti as a world-tube sliced by a family (St) of MOTS. If, 
in addition, the condition {8} is enforced, then IH inner 
BCs pi, HO, HI IH H are recovered in the stationary 



limit and one of our basic requirements is fulfilled. Even 
though the excision world-tube Ti is indeed expected to 
be a DH in realistic contexts, such a character is not 
actually enforced since the MTS condition 9^ < is 
not explicitly imposed. This is not a shortcoming of 
the approach. In fact, an (arbitrary) negative value for 
6^ k ' could be explicitly enforced as a Robin condition on 
/3 X -* 2 (cf. Eq. (16)inRef. [H): together with Eqs. ©, 
([6]) this would fix N, therefore providing an alternative 
manner of fixing the slicing. However it is known that 
the future evolution of a DH can cease "momentarily" to 
satisfy MTS and FOTH conditions, e.g. in the merging 
of two black holes once the common horizon has shown 
up. In this situation, insisting in the prescription of a 
negative probably leads to the ill-posedness of the 
whole coupled elliptic-hyperbolic system. For this rea- 
son, we rather adopt the methodological choice of only 
prescribing MOTS as inner BCs. Regarding a possible 
FOTH condition failure, and according with the charac- 
teristic analysis in pj], monitoring the sign of (/3 1 - — N) 
determines if inner BCs must or must not be provided 
for the radiative modes. This work represents an inter- 
mediate step in the ongoing program [6J addressing fully- 
constrained excised black hole numerical evolutions. 
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